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Abstract. It is shown for two electron atoms that ground-state wavefunctions 
of the form 

= </)(ri)</i(r 2 )(coshari -|-coshar 2 )(l-h 0.5ri2e“*"^i^) 
where rl and rj are the coordinates of two electrons and ri 2 = |ri — r 2 |, can 
be made highly accurate by optimizing a, b and tj>. This is done by solving a 
variationally derived equation for cji for a given a and b and finding a and b so 
that the expectation value of the Hamiltonian is minimum. For the set {a, b, 0} the 
values for various quantities obtained from the above wavefunction are compared 
with those given by 204-parameter wavefunction of Koga et al [ll]and are found 
to be matching quite accurately(within ppm) with them. 


1. Introduction 

The simplest of many-electron systems are those with two electrons like ions of the 
helium isoelectronicc series or a harmonium atom. Since the Schrodinger equation for 
these systems cannot be solved analytically, accurate solutions are obtained by using 
the variational principle for the energy. Thus an approximate parametrized form of 
the wavefunction is chosen incorporating the effect of electron-electron interactionon 
in the wavefunction and the expectation value of the Hamiltonian (atomic units are 
used) 

H = ——Vj — — V 2 + Vext{l^l) + Vextil^) H- (1) 

2 2 ri2 

is minimized with respect to the parameters in the wavefunction. This leads to an 
optimized approximate wavefunction and the corresponding energy. The various forms 
used for the wavefunctions are those by Hylleraas[l, 2, 3, 4, 5], Chandrasekhar [6] and 
Kinoshita[7, 8]. As an example of how these forms incorporate the effect of electron- 
electron interaction, we consider the Chandrasekhar wavefunction for two-electron 
atom given as 

^c(rT,r5) = (2) 

where a and b are the parameters to be optimized and Ctv is the normalization 
constant. The idea behind writing the wavefunction in the form above is as follows. In 
the ground state of a two electron atom, each electron moves essentially in a hydrogen- 
atom like Is orbital (of the form e““’’)but because of electron closer to nucleus, the 
effective nuclear charge is different for the outer electron. Hence the parameters 
for the orbitals of the two electrons are different. Symmetrization of the resulting 
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product wavefunction then leads to the form given in equation (2). Another way of 
obtaining accurate wavefunctions has been given by Hylleraas[l, 2, 3, 4, 5] where the 
wavefunction is expressed as 

N 

(rl, rl) = ^ (3) 

i=l 

where s = ri +r 2 , t = ri — r 2 , u = \ f{ —r^|, and C, and Ci are variational parameters. 
The most accurate calculation using this form of the wavefunction has been that of 
Frankowski and Pekeris[9], who used a 246-term wavefunction to get highly accurate 
energies for the helium isoelctronic series. Employing a 230-term wavefunction Freund 
et al. [10] got essentially the same results as those of [9]. Using another modification 
to Hylleraas type wavefunction, Koga et al [11] constructed a 204-term wavefunction 
that gives energy within 16 nano-Hartree of the energy given in [10]. In addition, 
usefulness of the work of [11] lies in the fact that they expressed the electronic density 
of the 204-parameter wavefunction in a simple analytical form with 43 parameters. 
Along similar lines, in 1994 Umrigar and Gonze[12] constructed very accurate densities 
for He isoelectronic series by using 491-term wave function having the same form of 
the basis set as used by Freund et al . The resulting densities have been used in the 
past to construct exchange-correlation potential[12] and to calculate[13] polarizability 
of two electron atoms and also to study derivative discontinuity of the exchange- 
correlation[14] of potential of density functional theory. Examples of references 
[12, 13, 14] are sufficient to show that the availability of accurate densities facilitate 
studies in fundamental aspects of density functional theory [15]. As such methods that 
can give accurate densities with relative ease are desirable to facilitates such studies. 
With this in mind Le Sech[16] introduced a semianalytical wavefunction for the ground 
state of two electron atoms and ions that gives energy within parts per million (ppm) 
of the answers of references [9, 10, 11]. 

2. Le Sech wavefunction 

The wavefunction proposed by Le Sech is motivated by the intuitive aspects of the 
Chandrashekhar wavefunction and correlation factor given by earlier studies [17]. The 
wavefunction is given as 

^L(bi,b2) = C]\[e~^^''^~^'~^\coshari coshar2)(l -I- 0.5ri2e“^’’^^) (4) 

where Cjv is the normalization constant and Z is the atomic number of the atomic 
system. The factor (coshari -|-coshar 2 ) is chosen so that cusp condition[18] is satisfied 
exactly by the wavefunction. The parameters a and b are obtained variationally by 
minimizing the expectation value of the Hamiltonian with respect to them. The 
energies so obtained for H', He, Li+, Be^+j are given in Table 1 along with the 
exact energies obtained in [11]. The differences between the two energies are also 
shown. We mention that the values of parameters a and b and the corresponding 
energies obtained by us are slightly different from those reported in ref [16]. It is clear 
from Table 1 that the wavefunction 412, of Equation (4) leads to accurate energies with 
the difference from those of [11] becoming smaller with increasing Z. The difference 
for H“ is 1899 parts per million and goes down to 104(ppm)for We note that the 
wavefunction of Equation (4) can be improved further with the inclusion of one more 
parameter and leads to a closer agreement with the exact energies. 


Highly accurate wavefunctions for two-electron systems using two parameteres 


3 


Table 1. Energies for He-like systems using approximate wavefunction of 
Equation (4) and of Equation (9). Their comparision with the energies given 

in [11] (given in the column under’Ifx) is also made in the table. The differences 
of each energy from that given in [11] is given in ppm under the energies. 


Atom 


■^K 



a 

b 

-Energy(a.u.) 

-Energy(a.u.) 

-Energy(a.u.) 

a 

b 

H- 

0.58 

0.06 

0.5267 

0.5277 

0.5271 

0.62 

0.06 

He 

0.72 

0.20 

(1899) 

2.9020 

2.9037 

(1137) 

2.9028 

0.93 

0.20 

Li+ 

0.87 

0.36 

(586) 

7.2778 

7.2799 

(310) 

7.2788 

1.19 

0.36 

Be2+ 

0.99 

0.52 

(289) 

13.6533 

13.6555 

(151) 

13.6544 

1.48 

0.54 

B^+ 

1.10 

0.67 

(161) 

22.0286 

22.0309 

(81) 

22.0297 

1.72 

0.70 




(104) 


(55) 




The question we now ask is how accurate are the densities given by '1'^? This is 
important from the point of view of having an accurate and easily accesible density 
of two electron systems if we were to employ them to investigate the foundational 
aspects of density functional theory of many electron systems in general and two elec¬ 
tron systems in particular. To this end we present in Table 2 the values of p(r = 0), 
^|r=o and in Table 3 and 4 various moments of the densities obtained from the Le 
Sech wavefunction. A comparision of the results is made with the accurate values of 
those obtained from the densities given in [11]. It is seen that the densities given by 
are smaller than the exact values of the densities at r = 0. Since the Kato-cusp 
condition is satisfied by both and the derivative |^|r=o are also small in 
comparision with the exact density. This implies that the Le Sech densities may be 
little more spread out than the exact densities. We also test the density in terms of 


Table 2. Densities and its derivative at r = 0 are calculated by using ^ml- 
For the density given by Koga et al (given in the column under i/jk) we have used 
the exact expression given in their paper[ll] to calculate p(r = 0) and —■^(^" = 0)- 
We see that Le Sech density and its derivatives are smaller than the density given 
in[ll] at r = 0 but the Modified Le Sech density is close to them. Numbers given 
are in atomic units. 


Atom 

II 

o 






^ ML 



ML 

H- 

0.310 

0.323 

0.329 

0.621 

0.649 

0.657 

He 

3.554 

3.621 

3.616 

14.214 

14.483 

14.462 

Li+ 

13.552 

13.704 

13.693 

81.312 

82.223 

82.185 

Be2+ 

34.180 

34.396 

34.384 

273.438 

275.167 

275.077 

B3+ 

69.164 

69.517 

69.501 

691.640 

695.165 

695.011 


it satisfying the ionization potential (IP) theorem in Kohn-Sham density functional 
theory. Furthermore, we also study the accuracy of the corresponding Kohn-Sham 
exchange-correlation potential. According to the IP theorem [19, 20] the highest oc- 
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Next in Table 3 and Table 4, we present various moments of the density obtained 
from d'L of Equation (4) and compare them with their exact counterparts. It is seen 
that (r“^) and (r~^) for are less than those obtained from[ll] the exact density. 
This is consistent with the density for T ^ being smaller than the exact one near the 
nucleus. For (r), (r^) and (r^) on the hand, it is observed that their values are larger 
than the exact moments for most of the cases (exception is for H'). This again shows 
that for 'I'l, the density is more spread out in general. 

Table 3. Comparision of moments calculated by using of Equation (4) and 
'^ML of Equation (9) with the moments given in [11] for two electron atoms. The 
latter are given in column under ^k- Numbers are given in atomic units. 


Atom 

(r-q 

(r q 



•^K 

ML 



ML 

H- 

2.145 

2.233 

2.229 

1.348 

1.366 

1.367 

He 

11.878 

12.035 

12.002 

3.363 

3.377 

3.376 

Li+ 

29.574 

29.855 

29.762 

5.362 

5.376 

5.375 

Be2+ 

55.301 

55.680 

55.485 

7.367 

7.376 

7.376 

B3+ 

88.900 

89.507 

89.146 

9.366 

9.375 

9.375 


Table 4. The caption is same as that of Table 3 


Atom 

(r> 

(r^> 

(r^) 


'i’L 

•Hk 

4' ML 

'^L 

4'k 

I'ml 

4'l 

'i'K 

4'ml 

H- 

5.331 

5.420 

5.422 

22.000 

23.830 

23.859 

125.030 

152.000 

152.366 

He 

1.863 

1.859 

1.859 

2.402 

2.387 

2.386 

3.996 

3.936 

3.934 

Li+ 

1.148 

1.146 

1.145 

0.898 

0.893 

0.892 

0.895 

0.883 

0.882 

Be^+ 

0.829 

0.829 

0.829 

0.465 

0.464 

0.464 

0.331 

0.328 

0.328 

B3+ 

0.650 

0.649 

0.649 

0.285 

0.284 

0.284 

0.157 

0.156 

0.157 


cupied orbital eigenenergy is equal to negative of the ionization potential of a many 
electron system. Accordingly for a Kohn-Sham system for a given ground-state den- 
sity of two electron system, its eigenvalue should be equal to / = — (E + ^) where E 
is the energy of the two-electron system. To test the accuracy of the Le Sech densities, 
we construct the corresponding Kohn-Sham system using the Zhao-Parr[21] method, 
which is described briefly in the following. 

For the orbitals of noninteracting system of electrons the kinetic energy Tg and 
ground state density p(r) are given by 

T’s = - ^V2|(/)i(r)), p(r) = y] (5) 

i i 

For a given ground state density, the Kohn-Sham system is formed by minimizing 
Tg + J Vext(r)p(r)dr with respect to ffs with the constraint. 

[p(r) - po(r)][p(rO - po(rO] ^ ^ 

|r - r'l 



( 6 ) 
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which is equivalent to p{r) — pg(r) = 0, and / |(^ipdr = 

- + Vext{r)(t)i{v) +v{r)(j)i{Y) 


1. This leads to the equation 
= Ci(fr{Y) (7) 


where 


7(r) = xj 


[p(r')-Po(r')] 


dr' 


( 8 ) 


where A and are the Lagrangian multipliers to enforce the two constraints above. 
Equation (7) is solved self-consistently using a Herman-Skillman[22] program modified 
suitably. Since at the solution point p(r) = po(r), the value of A should go to infinity 
such that v(r) remains finite. In this paper we have performed all the calculations for 
A = 30000. We point out that for two electron systems, the Kohn-Sham systems can 
also be found directly from the Laplacian of the density. We however, chose Zhao-Parr 
method because it is more general and more accurate. 


Table 5. Highest occupied orbital eigenvalue emax and ionization potential 
I = —{E H—for two electron systems. Numbers given are in atomic units 


Atom 


L 

'^ML 

Iexpt[23] 


^max 

I 

^max 

/ 


H- 

0.0613 

0.0267 

0.0276 

0.0271 

0.0277 

He 

0.8265 

0.9020 

0.9059 

0.9027 

0.9036 

Li+ 

2.5718 

2.7778 

2.7821 

2.7787 

2.7798 

Be^+ 

5.3581 

5.6533 

5.6593 

5.6543 

5.6557 

B3+ 

9.1406 

9.5286 

9.5349 

9.5296 

9.5320 


Shown in Table 5 are the values of —Cmax and / for for the ions studied in 
this paper. It is seen for dt l that the values of Cmax and / are comparable to each 
other. However, they are not close enough to satisfy the IP theorem. Furthermore, 
the values of / are close to the experimental ionization potential, which shows that 
the total energies E for d'i is accurate. On the other hand —Cmax is not accurate 
for This is particularity evident for the densities of the H~ ion and He atom 
where the difference between the two is 121 and 8, percent respectively. Since the 
eigenvalue Cmax is determined by the asymptotic decay of the electron density, it is 
evident that asymptotically, the density given by Le Sech wavefunction is not accu¬ 
rate. This is further evidenced by the corresponding Kohn-Sham exchange-correlation 
shown in figure 1 where we have plotted the exchange-correlation potential of He for 
both Le Sech density as well as the density given in[ll] and another accurate density 
of Umrigar and Gonze[12]. It is clear that Le Sech density leads to smaller magnitude 
of the exchange-correlation potential near the nucleus, which in turn gives densities 
which are more spread out. 

It is evident from the results presented so for that the Le Sech wavefunction 
needs to be improved if it were to give accurate densities. The rest of the paper is 
then devoted to modifying the Le Sech wavefunction so that it leads to densities that 
can be employed to study exact properties in density functional theory. 
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Figure 1. Comparision of the exchange-correlation potential calculated for 
given Le Sech density (Vxcl), modified Le Sech density (Vxcml) for 

the density given in[ll] (Vxcx) and for the accurate density of Umrigar and 

Gonze[12]{V xcug)- 



Exchange-correlation potential for He 



r(a.u.) 


3. Modified Le Sech wavefunction 

In this section we modify the Le Sech wavefunction so that the new wavefunction gives 
accurate energies as well as densities. We first point out how the Le Sech wavefunction 
can be made better. We feel that keeping the parameter Z (of the factor) fixed in 
the Le Sech wavefunction does not provide enough flexibility in the wavefunction and 
it is this component of the wavefunction that needs to be calculated more accurately. 
Hence the modified Le Sech (ML) wavefunction for the ground state of two electron 
atoms that we propose is 

'E'ML(hi,r2) = (/)(r-i)(()(r2)(coshar-i + coshar2)(l + 0.5ri2e“'"’i^) (9) 

where for a given set of a and b, the orbital (j)(r) is to be calculated self-consistently 
by solving the equation derived below which is obtained from the Euler equation 

\ IW)mL — E{'i>ML\'^ML)] = 0 ( 10 ) 

d(p[rj 

where {H)ml is the expectation value of the Hamiltonian of equation (2) with respect 
to the wavefunction '^’ml above. Equation (10) leads to following Schrddinger-like 
equation for (in the equation below ri 2 = |f — r^l, V is gradient operater in 
terms of r and V 2 in terms of r^). 

- ■ J f(rH2,ri2)|(()(r2)pVf(r,r2,ri2)dr2 
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1 


■ J [|?^'(r2)Pf(r, r2, ri2)V2f(r, I2,112) + |(/>(r2)pf(r, r2, ri2)V^f(r, I2, ri2) 


2A{r) 

+ |f(r,r 2 ,ri 2 )P^(r 2 )V 2 <?i(r 2 ) + 2f(r, 12 , ri 2 )(/)(r 2 )V 2 (/>(r 2 ) • V 2 f(r, r 2 , ri 2 )]dr 2 (/)(r) 


r A{r) 


l</'(r2)f(i',r2,ri2)p 


r2 


■dr2<^(r) 


+ 


1 f l<A(i'2)f(r,r2,ri2)p 


A{r) 


r-r2 


■dr2</>(r) = Ecj){r) 


where 

and 


/(r, r 2 , ri 2 ) = (coshor + coshar- 2 )(l + 0.5ri2e 

A{r) = J |?i>(r2)f(r,r2,ri2)pdr2 


( 11 ) 

( 12 ) 

(13) 


The equation above is solved self-consistently for a range of parameters a and b. The 
solution gives 4){r) and the corresponding energy E as the eigenvalue. The best wave- 
funcion corresponds to the values of a and b that give the minimum value of E. 


We wish to point out that after completing our work, we discovered that the 
above equation was solved[24] in 1937 to obtain ground state energy and polarizabil¬ 
ity of the He atom. However the scope of the present paper is much wider. In this 
paper we investigate the wavefunction not only for the energy but also for the density 
it gives rise to. We present the results obtained by solving Equation (11) for the two 
electron atomic systems studied so for and show the accuracy of the energy as well 
as the density that is achieved by self-consistently determined Tml- From the re¬ 
sults presented below, it becomes clear that Equation (11) provides a computationally 
straightforward method of obtaining highly accurate wavefunctions for two-electron 
systems. The method can also be applied with equal ease to other two-electrons sys¬ 
tems such as electrons in a confined environment or harmonic oscillator potential. 


Given in Table 1 are the results for the total energies of the modified Le Sech 
wavefunction Tml of Equation (9). As is evident, for the modified wavefunction, the 
total energies are better than those for the Le Sech functions. The modified wavefuc- 
tion reduces the error of energies from Le Sech wavefunction by roughly a factor of 2. 
This is because the wavefunctions has been made more flexible. In addition to giving 
improved energies, 'i>ML also gives much better densities than Tl. We present these 
results in the following. 

In Table 2, we also present the density and its derivative at r = 0 obtained from 
d^ML of Equation (9). We see that the values of modified Le Sech density and its 
derivative are very close the values obtained using exact expression given in [11]. In 
Table 3 and Table 4 we present the various moments of the density obtained by 'it ml- 
These moments are compared with the moments given in [11]. We get all the moments 
calculated by the modified Le Sech density close to the moments given in [11]. This 
implies that the modified Le Sech density is very close to the exact density. 

In Table 5, IP theorem is tested for both Le Sech and modified Le Sech density, 
where the e-max are calculated by solving Equation (7). We see that the modified Le 
Sech density satifies the IP theorem almost exactly. Finallly in Figure 1 we show 
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the exchange-correlation potential by using Le Sech density, modified Le Sech density, 
density given in[ll] and for the density obtained by Umrigar and Gonze[12]. We see 
that the exchange-correlation potential given by modified Le Sech densities is very 
close to the exact potential. 

4. Conclusion 

Le Sech wavefunction[16], which gives very accurate energies for He-like systems, has 
been modified so that both the energies and densities obtained from the modified 
wavefunction are highly accurate. The densities so obtained can be used with 
confidence to perform fundamental density functional theory investigations. Along 
this line of investigation we are now employing the method proposed to study adiabatic 
connection in density functional theory. 
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